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ABSTRACT 



Context. Methanol masers at 6.7 GHz are recognised markers of high-mass star formation regions. The study of their distribution 
in the Galaxy gives important insights into the star formation activity of the Milky Way. We present a statistical analysis on 
the General Catalogue of 6.7 GHz methanol masers in the Galaxy with the aim of extracting global properties of the masers. 
Aims. We provide constraints on the luminosity function of 6.7 GHz methanol masers and on their total number in the Galaxy. 
Methods. We model the spatial distribution of the masers in the Milky Way by using their distribution in galactocentric 
distance which is unambiguous once a rotation curve for the Galaxy is assumed. This is the starting point for determining the 
luminosity function of the masers. 

Results. The luminosity function of 6.7 GHz methanol masers is modelled as a power-law with sharp cutoffs and having an 
index lying between —1.5 and —2. We also predict the number of detections of methanol masers assuming different sensitivity 
limits in the observations. 
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1. Introduction 



Since its discovery by iBatrla et ail (|l987ri and iMentenl 
( 19911 ) the bright maser emission from methanol (at 12.2 
and 6.7 GHz respectively) has become a reliable tool for 
detecting and studying regions where (massive) stars form 
and are in the ir very early stages of evolution (see e.g. 
Ellingsen 20061 ). Methanol masers were divided into two 
empirical classes, I and II (|Mentenlll99lf) . Class II masers 
are detected close to strong Infrared sources (as e.g. 
Ultra Compact (UC) Hn regions), while class I are ob- 
served offset from these objects, in the shock regions of 
their outflows. Theoretical modelling was able to identify 
class I masers as collisionally pumpe d masers, while class 



II as radiatively pump ed ones (e.g. iSobolev et al.1 11997 



Cragg et all [200lL [20051. 



Observations of methanol maser sites at other wave- 
lengths support the paradigm of methanol masers as 
one of the first signposts of massive star formation (e.g. 
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Walsh et al.lll999l : iGoedhart et~aT1 120021 : IPestalozzi eTal 
20023). Being very bright, methanol masers are ideal 
for high spatial resolution observations using interferom- 
eters where very detailed positioning as well as map- 
ping of the finest spatial and dynamical maser features 
gives importa nt insights into the natur e of some known 



sourc es (e.g. 



Minier et al. 2001, 2003: Pestalozzi et al 



20041 ). Nevertheless, the question of whether these masers 
trace discs or outflo ws in young pr otostars still remains 
open to debate (e.g. DeBuizerll2003l ). 

Their association with youthful and massive star birth 
opens up the possibility of using methanol masers as 
a new and reliable tracer of that rapid stage of evolu- 
tion prior to the development of UC H11 regions. For 
more than a decade now, searches for new methanol 
maser sour ces have been u n derta k en by a numbe r of au - 
thors, e.g. iMacLeod etail (Il992h: ISchutte et all (Il993l): 



Caswe 1 et all (Il995lh lEllingser d 19961 ): ISzvm czak et al . 
(|2000aL l20o3 ) : IPestalozzi et al.l (|2002bh . All together, this 
brought the number o f known methanol masers to 519 
(jPestalozzi et al.ll2005l ). This represents a statistically sig- 
nificant sample which, for this paper, motivates us to 
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study their spatial distribution throughout the Galaxy 
with the aim of determining their luminosity function. 

To date, most statistical work has focused on find- 
ing correlations between the physical characteristics of 
the maser and the associated IRAS source. In general, 
6.7 GHz methanol masers seem to be more efficiently de- 
tected towards bright IRAS sources, having fen > 100 Jy 
(|van der Walt et alJll996t ISzvmczak fe Kusll2000bl ). with 
a clearly hig her detection rate in the inner than in the 
outer Galaxy ( Szvmczak fe Kusll2000bl ). Early attempts to 
find a correlation between the associated IRAS flux den- 
sit y with the maser flux den sity gave no positive results. 
In van der Walt et all (|l995h the authors are left with the 
fact that the maser flux densities are smaller than the 
100/im ones, suggesting that the masers could be pumped 
by the lOO^m photons. The lack of correlation between 
maser flux density and IR flux density is not explained in 
that work. IRAS sources associated to 6.7 GHz methanol 
masers seem to be concentrated in a small region of the 
[25-12]-[60-25] colour-colour diagram, i ndicating that the 
maser arises in specific environments ( Szvmczak fc Kusl 
20001)). Based on the same conclusions, IXu et al.l ((2003) 
claim that the methanol maser phase must be short and 
occurs in the early stages of star formation, where the 
IR luminosity is high enough that the IR radiation it- 
self might be responsible for the ma ser pumping, in ac- 
cordance with the oretical modelling (jSobolev et al.lll997 : 



Gragg et all l2005h . One should note, however, that the 
association of methanol masers with IRAS sources is 
not a reliable one, a s several studie s at high resolution 



have confirmed (e.g. Ellingsen 19961 : Minier et al. 2002t 



Pestalozzi et al. 2002q ): methanol masers are often seen 



offset from all IRAS or centimetre continuum sources. 

The range of velocities and the lincwidths of methanol 
masers have also been used to try to deduce so me intrinsic 
charac teristics for these objects. For instance, Slvsh et "all 
1999) found that methanol masers show a clear velocity 



difference relative to the velocity of the parent cloud (a 
comparison made with the velocity of CS). This could sug- 
gest an origin for methanol mas e rs in di scs seen edge-on. 
Alternatively, ISzvmczak fc Kusl (|2000bl) suggested that 
the velocity dispersion seen in maser spectra could be used 
as an evolutionary tool, arguing that the early stages of 
the formation of an UC Hll region would be traced by 
narrow line masers because of the low number of mas- 
ing clouds and the lower velocity dispersion within them. 
Once the UCHII region begins to disperse the surround- 
ing circumstellar matter, larger linewidths are expected. 
All such hypotheses which appeal to linewidths or the ra- 
dial velocities of spectral features are subject to the caveat 
that spectral feature blending could falsify the measure- 
ment of the lincwidth resulting in, for example, a false 
evolutionary sequence. 

Comparisons hav e also been made between masers at 
different frequencies. Slvsh et al. ( 1999h find that in those 
sources where 6.7 and 44 GHz masers coexist, the ratio 
of their flux density can be used to discriminate class 
I from class II masers. On the other hand, no correla- 



tion was found when comparing 6.7 with 12.2 GHz coun- 
terparts, except for the clear fact that 6.7 GHz masers 
are always found to be stronger in flux and in bright- 
ness temperature when ass uming a constant source size 
()Malvshev fc Sobolevlliool . 



More recent studies have concentrated their efforts 
in c haracterising the sources hosting methanol masers 
(e.g. lHill et all2005t|Purcell et al.l2006fh Methanol masers 
seem to be always associated with 1.2 millimetre dust 
continuum emission as well as CH3CN emission. Maser- 
hosting sources are in general more massive and hotter 
than sources showing 1.2 mm continuum emission but no 
maser. Also, sources hosting methanol maser show hot 
core characteristics and most of them pass the Lumsden 
MSX criterion for massive protostars in the Galactic 
plane (jLumsden et al.l 120021 ). Finally, methanol masers 
have been detected toward dark mid-IR clouds, clearly 
indicating that star formation is ongoing in those regions. 

The luminosity function of 6.7 GHz methanol masers 
has been studied using different approaches. The main 
problem is in determining the heliocentric distance for ev- 
ery source. Because of the high corre spondence of 6 . 7 GHz 



methanol masers with OH masers, ICaswell et al.l (J1995) 



assume that the luminosity function of 6.7 GHz methanol 
masers should be simila r to the one for OH masers (see 
ICaswell fc Havneslll987l) from which they conclude that 
there must be some 500 meth a nol m aser sources in the 
Galaxy. In van der Walt et al. (|l996l ). the ambiguity of 



kinematic distances is solved in a probabilistic way. When 
the decision was not clearly made on the basis of the 
total luminosity of the IRAS source hosting the maser, 
the source was assigned a probability for it to lie at the 
near heliocentric distance. Repeating this assignment in 
a random way (and for different probabilities) and av- 
eraging over the results the authors obtained a spatial 
distribution from which the luminosity function was es- 
timated. For their sample of about 240 sources, the au- 
thors fit a power-law luminosity function with an index of 
rj — 2 (the fitting was done on the distribution of sources 
per unit luminosity interval). This index is expected to 
flatten in the low-luminosity end of the distribution as 
otherwise the total number o f methanol mas ers in the 
Galaxy would be far too high. lElhngsen (|l996f) refrained 
from any luminosity function estimate, invoking the prob- 
lem of determining the distance to every source. More re- 
ce nt studies of large sample s of 6.7 GHz methanol m asers 
bv ISzvmczak fc Kusl (|2000bh : lsivmczak et all (|2002l ) have 
not directly addressed the study of the luminosity function 
of methanol masers. 

The most recent statistical study of methanol masers 
in the Galaxy has been presented bv lvan der Walt (2005) 
in which estimates for the lifetime of the methanol maser 
phase and the total number of methanol masers are made. 
The model uses a combination of the initial mass func- 
tion, (local) star formation rate, a synthetic distribution 
of these sources in the Galaxy and a constant detection 
limit of 1.5 Jy in some template regions to conclude that 
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methanol masers should last between 2.5 — 4.5 x 10 4 years 
and that there should be around 1100 in the Milky Way. 

The aim of this paper is to present a statistical anal- 
ysis of the population of 6.7 GHz methanol masers in the 
Milky Way. We aim to constrain the luminosity func- 
tion, and begin by modelling the spatial distribution of 
methanol masers in our G alaxy. This approach is diff erent 
from the one presen ted in Ivan der Walt et aD (|l996l ) and 
van der Waltl (|2005i ). where the authors apply a random 



choice algorithm to solve the heliocentric distance ambigu- 
ity for sources on galactic orbits internal to the Sun. Our 
starting point is the distribution of methanol masers in 
(kinematic) galactoccntric distance, which is unique once 
a rotation curve for the Galaxy is assumed. We then as- 
sume that that distribution is the azimuthally averaged 
surface density of sources in the Galaxy. In this way we 
are not compelled to solve the distance ambiguity for ev- 
ery source lying on galactic orbits internal to the Sun, as 
the decision is made by the surface density itself. 

Nevert heless, the es t imate of the total number of 
sources in van der Ward (2005) will be used as check for 
consistency in this work. Furthermore, the present study 
repre sents a study of pri nciples with fewer assumptions 
than Ivan der Waltl (|2005h . Finally, the physical meaning 
of the luminosity function of the maser emission is not 
addressed in this paper. 

2. Modelling the maser distribution and 
luminosity function 

The determination of the luminosity function for astro- 
nomical sources requires an accurate knowledge of their 
distance from the observer. In the case of galactic sources, 
kinematic distances from spectral line observations are be- 
coming more accurate thanks to improved modellin g of the 
Galactic rotation curve (e.g. iBrand fe Bli"t3ll993l as well 
as iRusseill 120031) . The difficulty resides in the discrimina- 



tion between near- and far- heliocentric distances for all 
sources on orbits internal to the Solar orbit around the 
Galactic Centre. 

Beside the probabil istic approach presented in 



van der Walt et al 



( 1996L another way to solve the he- 
liocentric distance ambiguity is the observation of Hi self- 
absorption. The method relies on the fact that molecu- 
lar clouds have an outer layer of neutral hydrogen which 
can absorb the ubiquitous background galactic Hi emis- 
sion. Sources lying on the near side of the galactic centre 
would have more Hi emission to absorb than sources on 
the f ar side, making it straightforward to solve the ambigu- 
ity ( Jackson et al.112002 ; Liszt et aTlll981 ) . As the previous 



one, this method still relies on the calculation of kinematic 
distances, which requires a goo d knowledge of the ro tation 
curve of the Galaxy. Recently, iBusfield et al. (2006) used 
this method to successfully determine the distances of a 
number of massive young stellar object candidates in the 
fourth quadrant of our Galaxy. The authors succeeded to 
resolve the distance ambiguity for 80% of their targets, 
the main limitations being the coarse spatial resolution of 



the archival Hi data as well as the inability to apply the 
method for sources at galactic longitudes £ < | 15° |. 

The main idea of the present study is to avoid solving 
the distance ambiguity for every single source and adopt 
a global approach instead. This approach starts with a 
model of the spatial distribution of methanol masers in 
the Galaxy. We then model the observability of sources 
with different luminosities at different detection limits and 
estimate their real total number. 

The sample of methanol masers used for the present 
study is the General Catalogu e of 6.7 GHz Methano l 



Masers (GCMM) published by IPestalozzi et all (|2005l ). 



available on the CDS websit^0. This sample is far from 
complete, but allows a good study of principles and a first 
test for the formalism we adopt. The outcomes of this 
study will be compared to the ones from the Methanol 
MultiBeam Survey (MMB Survey % The MMB Survey is 
expected to probe the whole Galaxy at a depth of 0.1 Jy 
at lcr, and will probably reveal the complete population 
of 6.7 GHz methanol masers in the Galaxy. 



2.1. Spatial distribution of the masers 

For a given observational sensitivity, the observed spatial 
distribution of a certain type of object in the Galaxy is 
strongly dependent on the luminosity function of those ob- 
jects. If, for instance, the luminosity function is strongly 
peaked towards low luminosities, all searches for those ob- 
jects at a certain detection limit will select the majority of 
sources in the solar neighbourhood and miss those lying 
further away. This will give an observed spatial distribu- 
tion of sources which will probably not reflect the real 
distribution of sources. If on the other hand the spatial 
distribution is known to a high degree of confidence it is 
then possible to estimate both the luminosity function as 
well as the real total number of those sources. 

A general view of the methanol masers known to date 
superimposed on the CO emission in the galactic plane 
both in space and LOS velocity is presented in Fig. [TJ 
The masers clearly follow both the structural and dynam- 
ical features of the CO gas in the Milky Way. There is a 
clear concentration of sources at longitudes close to ±50° 
which also follows the dynamical signature of the molec- 
ular ring. Also, the mase rs seem to be concen trated in 
the galactic plane (see also IPestalozzi et"ld 1 l2005l Finally, 
there are some masers that seem to mark the nuclear ring 
(I ~ 0°), and hence lie within 1 kpc from the Galactic 
Centre. This is not the first time that masers are detected 
within 1 kpc of the Galactic Centre (see e.g. the detection 
of OH masers, Caswell fe Havnesl 1983 ). The discussion 
on this point goes beyond the scope of this paper, but it 



1 http:/ /vizier, cfa. harvard. edu/viz-bin/VizieR?- 
source= J /A+ A /432 /737/ 

2 The MMB Survey aims at surveying a strip of b —\ 2° | 
across the galactic plane and at all longitudes, searching for 
6.7 GHz methanol masers. It has started in January 2006. See 



also http : / / www. j b . man . ac . uk/research/methanol / 
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Fig. 1. Distribution of methanol masers in Galaxy, superposed on the CO contours from lDame et al. I (Il987h . in space 



(top) and LOS velocity (bottom). The methanol masers seem to accurately follow the overall structure of the Galaxy, 
both in space and LOS velocity. Particularly visible in the bottom panel is the fact that methanol masers are tracing 
the spiral arms (150° > I > 80° and -40° > I > -90°) and the high rotational velocity of the nuclear ring (JwO"). 



is nevertheless interesting to note that, being methanol 
masers exclusively associated with star formation activ- 
ity, these detections are strong indications of recent star 
formation very close to the Galactic Centre. 

The starting point in our analysis is the distribution of 
masers against their galactocentric distance, as this does 
not suffer from any distance ambiguity. The distance of ev- 
ery source from the galactic centre is uniquely determined 
by its galactic longitude and LOS velocity, once a rotation 
curve fo r the Galaxy is assu med (we adopt the rotation 
curve of lBrand fc Blit zHl99#l . Figure [2] shows the surface 
density of methanol masers in the Galaxy (main panel) 
and the surface densities of methanol masers and molec- 
ular gas (inset, histogram and dashed line respectively) 
normalised to the integral under the curves. The ratios 
between the values at 5 kpc (peak, or molecular ring) and 
the values around the Sun (ss 8-9 kpc) are 6:1 and 5:1 
for masers and gas respectively. Knowing that methanol 
masers are associated exclusively with star formation re- 
gions, this fact could support the idea that star forma- 
tion is more efficient in the molecular ring than in the 
outer Galaxy. We do not draw this conclusion, as this ef- 
fect could be due to a lower maser detection rate com- 
ing from an uneven observational coverage of the outer 
Galaxy. Another interesting empirical fact is given by 
the ratio of the gas surface density to the maser surface 
density. This ratio is an estimate of solar masses of gas 
per methanol masers and it has values between 1.5 x 10 6 
and 9.7 x 10 7 for peak and solar neighbourhood, respec- 



tively. These values are comparable to the gas content in 



3 Note that for the calculation of the galactocentric distance 
we use the line-of-sight v elocity of the bri g htest spectral com- 
ponent. As mentioned in iPestalozzi et ail (|2005l ). the error of 
that calculation is about 1 kpc when taking into account the 
inaccuracies in the position and in LOS velocity. 



large molecular clouds (w 10 7 M Q , see e.g. iKnapen et al 
1993), which means that in the ring we are expecting 
sa 10 methanol masers per large molecular cloud, while 
in the outer Galaxy this number drops to approximately 
one maser every 10 clouds. In the same line of thought we 
can do an order of magnitude comparison of these numbers 
with th e total conte nt of gas in the Milky Way (f» 10 9 M G 



see e.g. |Pamdll993l) and obtain an empirical estimate of 
the total number of methanol masers in our Galaxy, which 
should be of the order of 10 3 sources. 

We consider the histogram in the main panel of Fig. [2] 
to be a sufficiently reliable signature for the shape of the 
spatial distribution of methanol masers in the Galaxy. We 
support our assumption with the fact that the masers 
seem to accurately follow the distribution of CO in the 
Milky Way both in space and LOS velocity (Fig. [T]). 
We make two assumptions: that the surface density of 
methanol masers in the Galaxy, F(x,y), can be projected 
onto the Galactic plane (| b \— 0°), and that the masers 
are distributed axisymmetrically, i.e. the surface density 
distribution depends only on the galactocentric radius, R. 
The total number of sources in an annulus of thickness 
cLR is then F(R) 2ir RdR = H(R)dR, where H(R) is the 
function fitted to the histogram in Fig. [5] We support 
our two assumptions with the following arguments: a) the 
distribution of methanol masers in galactic latitude has a 
very small FWHM (about 1.0°, see IPestalozzi et alJliooBh : 
and b) axisymmetry is assumed by noting that the clear 
peak in the histogram in Fig. [5] is at a galactocentric dis- 
tance where the spiral pattern of the Galaxy is n ot yet 
clearly visible (molecular ring, see e.g. Russeil 20031 ). The 
latter argument is also supported in iFuxl (|1999T ). where it 
is shown that the existence of a bar in our Galaxy does 



Pestalozzi et al.: Maser catalogue, the statistics 



5 



a> 
-a 

CD 



Galactocentric distance [kpc] 



Galactocentric distance [kpc] 

Fig. 2. Main panel: Spatial density of methanol masers in 
the Galaxy as a function of galactocentric distance. Inset: 
normalised profile of the surface density of the masers (his- 
togram) compared with the surface de nsity profile of the 
H 2 gas (dashed line, from iBlit d Il996h . See Table ED for 
the results of the fits to the histogram in the main panel 
(excluding the first bin). 



not significantly affect the dynamics outside the molecular 
ring. The total number of observed sources in the Galaxy 
Ntot-obs is given by the integral from to R ma x of the 
function H{R). The Gaussian profile used for the fitting 
is parametrized by the integral under the curve N, the 
mean R pea k and width a. The results of the fit are shown 
in Table [[J 



2.2. Luminosity distribution of methanol masers 
2.2.1. Step 1: equal luminosity for all masers 

We use the assumption of the masers having all the same 
luminosity in order to better understand the basic trends 
and behaviour of the model within a simple framework. 

If we assume that all methanol masers have the same 
luminosity, the sensitivity limit of the observations trans- 
lates into a maximum heliocentric distance r ma x-obs, 
within which masers will be detected (dashed arcs a, b 
and c in Fig. [3]). We fit H(R) to the histo gram of Fig. [2] 
when applying different r max - bsi i.e. with a spatial dis- 
tribution which is truncated due to sensitivity. Recalling 
our basic assumption that the histogram in Fig. [5] is a re- 
liable signature for the shape of the galactic distribution 
of masers, we fix R pea k and a of the function F(R) found 
from the original fit to the histogram (Table [T]) and fit 
only for N. This will be the total number of sources in the 
Galaxy. We have then: 



where R m in/max arc defined by the heliocentric radius 
fmax-obs and the function max {R) ensures that masers 
are only counted if they lie within r max -obs- Note that 
if r-max-obs > Rq (i.e. the depth of the observations 
reaches beyond the Galactic Centre) then 9 max = ir for 
< R < r max -obs — Ro (arc a in Fig. [3|). The results of 
the fits are summarised in Tableland in Fig.[U We choose 
7 kpc as minimum value for r max -obs, as shorter distances 
would be strongly inconsistent with the observed distribu- 
tion of masers in the Galaxy. This is visible from Fig. [2l 
where the histogram shows that there are sources close to 
the galactic centre, i.e. at 8.5 kpc from the Sun. The fits 
are indicative of how many sources there should be in the 
Galaxy if the sample was severely limited by sensitivity. 



maser distribution 




N, nt - 



tot — obs 



+e max (i?') r R ma 



n(R') JR 



F{6, R') R' dR'd9 (1) 



Fig. 3. Geometry of the modelled methanol masers in the 
Galaxy. R pea k is the radius at which the maser distribution 
peaks, Rq the radius of the solar orbit around the galactic 
centre, 9 the central angle, r the heliocentric distance and 
I the galactic longitude. The dashed arcs a, b and c are 
equidistant curves from the Sun, described in the text. 



Increasing r max ~obs to large values (e.g. 40 kpc) is 
equivalent to stating that the sample in GCMM does not 
depend on the luminosity distribution of the masers, and 
hence that all masers in the Milky Way are detected. 
There are several reasons to believe that GCMM does not 
contain all methanol masers sources in the Galaxy: not 
only is there a bias due to the unknown luminosity func- 
tion but also the catalogue is biased due to non-uniform 
coverage of the galactic plane. The present study is meant 
to make predictions on the basis of the data available. 

Note that, by introducing a sensitivity limit, i.e. 
fmax-obs, we define a minimum luminosity for all de- 
tected masers. For a detection limit of 1 Jy over 0.2 kms -1 
or 5.56 kHz at t he rest frequency of 6.668519 GHz 
(jMuller et al.ll2004T ). we get minimum luminosities of the 
masers of 2.85 x 10" 6 , 1.91 x 10~ 7 and 7.33 x 10~ 8 L Q 
for r max - bs equal to 40.5, 10.5 and 6.5 kpc, respectively. 

Comparing these results with recent estimates of the 
total number of methanol masers in the Galaxy (as e.g. 
in Sect. 12.11 and van der WaltJ 2005f ). and retaining our 
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Table 1. Results of the best fit to the histogram in Fig. [21 
The values of r max _ f, s indicate a truncation of the func- 
tion F due to sensitivity. A large value of r max _ f >s (as 
e.g. 40kpc) is equivalent to no truncation (see text for an 
explanation). The column N lists the real total number of 
sources. The mean and width for the Gaussian distribu- 
tion are 5.08 and 1.42 kpc respectively. 
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Fig. 4. Dependence of the real total number of sources 
N inferred by fitting a distribution F truncated at differ- 
ent values of r ma x-obs to the histogram in Fig. [2] using 
the assumption of equal luminosity for all masers. The 
range of r m ax-obs defined by the vertical da shed lines cor- 



respo nds to the range for N estimated in Ivan der Walt 
2005. The horizontal dashed line indicates the total num- 
ber of sources in GCMM, 519. 



assumption that all masers have the same luminosity, we 
conclude that GCMM contains sources up to the galactic 
centre and just beyond (see dashed vertical lines in Fig.0]). 
The conclusions we can draw at this point of the paper 

are: 

— from fitting the histogram in Fig. [5] with truncated 
models of the spatial distribution (eq. [l} we see that 
GCMM has a considerable number of sources missing 
(up to a factor of about 2-3) ; 

— this suggests that the luminosity function of 6.7 GHz 
methanol masers is not a delta function (i.e. all masers 
do not have the same luminosity) and it is most likely 
dominated by intrinsically weak emitters, potentially 
young high-mass star formation regions; 

— the large number of missing sources is mostly due to 
the shallow sensitivity limit of most of the observations 
which contribute to the GCMM (see Pestalozzi et al 



2005). Only a very small part of the Galaxy was ob- 
served with a sensitivity limit better than 2 Jy. 



2.2.2. Step 2: introducing a distribution of luminosities 

The general distribution of luminosities of methanol 
masers is the link between the real spatial distribution and 
the observations. The observability of a source depends on 
a combination of its luminosity and the sensitivity of the 
survey used to detect it. High-sensitivity observations (e.g. 

0. Uy at la) can detect low-luminosity sources far away 
from the Sun, reaching a high level of completeness. Deep 
surveys guarantee that one probes most, if not all, of the 
Galaxy: this is the aim of the MMB Survey. 

By discarding the assumption of equal luminosity for 
all masers and introducing a luminosity distribution G(L), 
r m ax~obs no longer defines any sharp spatial cutoff point. 
The physical depth of the observations will scale as y/L/&, 
where L is luminosity and $ is the detection limit in the 
observations. It is therefore not possible to define an ab- 
solute value for the largest distance probed by the ob- 
servations, as this will vary within the same combination 
of luminosity function and sensitivity of the observations. 
For every sensitivity limit <1>, the scaling of the spatial 
cutoff becomes smaller increasingly quickly toward low lu- 
minosities. This means that if low luminosity sources are 
concentrated around some particular heliocentric distance 
r, these will be missing from the counts unless the sensitiv- 
ity is set accordingly. An order of magnitude in sensitivity 
implies a factor of y/lO in r max -obs- Assuming masers lie 
in a ring-like structure around the galactic centre (with 
some known radial profile as e.g. a Gaussian), such a dif- 
ference in r ma x—obs could mean either including or not a 
factor w2 sources in the final counts, as r max - b s could 
either include only the near peak or both the near and the 
far peaks of the distribution. 

To better understand the influence of the luminosity 
function in the present study we follow the following pro- 
cedure: 

1. Prescribe a functional form for the luminosity func- 
tion G(L). We require that this is normalised so that 
/ S{\.tl(L)) (ln(L)) = 1, where S(ln{L)) is the distribu- 
tion of luminosities sampled in equal logarithmic inter- 
vals. From this function, we construct a discrete set of 
probabilities Pj , centered on luminosities Lj , such that 

Zj Pj(Lj) = 1. 

2. Set a sensitivity limit $ for the observations. For every 
luminosity Lj , $ sets a maximum heliocentric distance 
fmax-obs(Lj, $) to which the observations will detect 
sources of that particular luminosity; 

3. For each luminosity bin we then: 

— Multiply the spatial distribution F(R) (without 
any distance cutoff) with P(Lj), i.e. the probabil- 
ity which corresponds to that luminosity. In this 
way we are left with the distribution of sources of 
luminosity Lj in the Galaxy; 
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— Apply the distance cutoff r max - b s (Lj, $) and 
count the number of sources in each galactocen- 
tric distance bin (this is equivalent to integrating 
over 6 and R as eq. [T]). In this way we count how 
many sources of the luminosity Lj we are able to 
detect with the sensitivity <I>. 

The total number of observed sources N to t-obs is given 
by the sum of all (detected) sources at all luminosities: 



N, n ,- 



tot — obs 



x Fi-t 



RdRdO (2) 



where _F ) _ trunc (i?) is the spatial distribution of sources 
within r max ^ b s (Lj, <£) and the integration limits of both 
integrals are the same as in equation [T] 

Once the shapes of the spatial distribution F(R) and 
of the luminosity function G(L) are chosen, the fit to the 
histogram in Fig. [2] has the following free parameters: the 
total number of sources in the Galaxy TV, the parameters 
defining G(L) and the sensitivity Note that we assume 
that the mean and width of the spatial distribution are 
fixed. 

To keep the number of free parameters to a minimum, 
we choose the luminosity function to be of the simplest 
kind, a power-law with sharp cutoffs L m i n and L max : 



G(L) =AL a 



(3) 



where A is determined by the normalisation condition, and 
the only free parameter defining G{L) is its power a. We 
take the range of luminosities f rom the literature. L m i n = 
1CT 8 L and L max = 10~ 3 Lq (jWalsh et al.lll997n . as our 
best estimate. 

2.2.3. Results 

From eq.[2]it is not possible to determine TV and a simulta- 
neously. Nevertheless, this method is still able to provide 
important insights and predictions on both the total num- 
ber of sources in Galaxy (TV) as well as on the slope a of 
the maser luminosity function. 

We can attempt to answer two different questions: 

— Knowing the mean detection limit in GCMM and 
how many sources GCMM contains, how many sources 
should there be in total in the Galaxy and how does 
this number depend on the slope of the luminosity 
function? 

— What fraction of TV do we expect to detect depending 
on the detection limit and the slope of the luminosity 
function? 

The answer to the first of these question is shown in 
Fig. [5] The curves in the graph were obtained by fitting 
equation [2] to the histogram in Fig. [21 excluding the inte- 
gral over R. For example, if we assume the mean detection 
limit of GCMM to be 1 Jy, TV would lie between » 520 and 
~ 5000, depending on the power of the luminosity func- 
tion. If we are to constrain our range for TV by taking 




Detection limit * [Jy] 

Fig. 5. TV versus detection limit at several values of a, 
the slope of the luminosity function (dashed -1; triangles 
-1.2; stars -1.5; dash-dot -2.0; circles -3.5). The spatial 
distribution profile used here is a Gaussian. The detection 
limit is defined as <& spread over 0.2 kms -1 . TV is the result 
of fitting eq. [2] to the histogram in Fig. [21 This explains 
the convergence to TV = 519 at low detection limits. 



CD 
W 
Q 




Detection limit <t> [Jy] 

Fig. 6. Ntot-obs as percentage of TV versus sensitivity limit 
at different values of a, the slope of the luminosity func- 
tion. The curves are labeled as in Fig. Note that an 
order of magnitude in sensitivity implies a factor of y/10 
in Tmax-obs- The spatial distribution profile used here is a 
Gaussian. The detection limit is defined as <& spread over 
0.2 kms^ 1 . 



into a ccount the predictio n made in Sect. 12.11 (TV « 1000) 
and of Ivan der Waltll2005l (TV w 800 - 1200), we can state 
that the luminosity function has a slope between ps — 1.2 
and « —1.5. At the detection limit of the MMB Survey 
(0.5 Jy) the total number of sources is expected to lie be- 
tween 519 and more than 2500, depending on the slope of 
the luminosity function. 
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The answer to the second question is shown in Fig. [5] 
In this approach the absolute number of sources N and 
Ntot-obs ar e irrelevant, and only the shapes of the spa- 
tial and luminosity distributions are important. From the 
graph we can say that by taking the mean sensitivity of 
GCMM to be 1 Jy, the total number of sources in that 
catalogue (519) represents between 10 and 80% of the to- 
tal depending on the chosen luminosity function. Again, if 
previous estimates of the total number of methanol masers 
in the Galaxy are correct, we can state that GCMM con- 
tains 30 — 50% of N, and this would constrain the slope 
of the luminosity function a to be > —1.5. The detection 
limit of the MMB Survey (0.5 Jy) will produce a catalogue 
containing between 15 and 85% of the real total number 
of sources in the Galaxy, depending on the slope of the 
luminosity function. 

Fig. [6] indicates that deeper observations have the ef- 
fect of reducing the dependence of the fraction N to t-obs/N 
from the luminosity function. A detection limit of O.Uy 
would allow to detect 60 — 95% of the total number of 
sources, which is a more favourable range than the one 
provided by a mean detection limit of 1 Jy. Assuming a 
power-law with slope —1.5, a 0.1 Jy detection limit would 
yield the detection of 85% of the total number of existing 
methanol masers in the Milky Way. 

3. Can we constrain the luminosity function from 
GCMM? 

Fig. [7] shows the distribution of luminosities of the masers 
in GCMM. Dots and diamonds represent the two extreme 
cases where all ambiguous sources have been assumed to 
lie at the near and far heliocentric distance, respectively. 
Since the sources suffering from distance ambiguity are 
about 90% of the catalogue, the two curves in FigUJshow 
a relative shift of about one order of magnitude in lumi- 
nosity. 

As GCMM is a highly heterogeneous sample (different 
sensitivities, non-uniform area coverage), we cannot con- 
clude that the plots in FigLHare indicative of the shape of 
the real luminosity function of the masers. Nevertheless, 
if GCMM does not contain most of the faint methanol 
masers in the Galaxy (< 1 — 3 Jy), we expect the counts of 
a large scale deep search of methanol masers in the Milky 
Way (as e.g. the MMB Survey) to populate the lower part 
of the luminosity range in Fig. [7] This is due to the fact 
that a 1 Jy source is intrinsically fainter than the peak lu- 
minosity in those distribution, as calculated in Sect. I2.2T1 
The detection of many sources at the 1 Jy level or be- 
low will change the shape of the luminosity distribution 
in FigUJ closer to one of a power-law. 

Maintaining the assumption of the luminosity function 
of methanol masers to be a single power-law between sharp 
cutoffs and observing it with a sensitivity of e.g. 2 Jy we 
obtain best similarity to the data as shown in FigUJ The 
synthetic data points (stars) seem closest to the data curve 
obtained when putting all ambiguous sources at the near 
heliocentric distance (dots). This strongly suggests that 



O 10 

O 
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x 



x 



Luminosity [Solar luminosities] 

Fig. 7. Distribution of the luminosities of all masers in 
GCMM, assuming, for all ambiguous sources, either the 
near (dots) or the far (diamonds) heliocentric distance. 
The luminosities were calculated by assuming that the 
peak flux is spread over 0.2 kms -1 . The difference in the 
peak position is of about 1 order of magnitude. Stars show 
the results obtained by the observation at 2 Jy sensitiv- 
ity of the modelled maser population in the Galaxy. The 
model is defined by a total number of sources equal to 
5000, a spatial distribution as in Fig. [2] and a luminos- 
ity function expressed as a single power-law of index -1.7 
between 10" 8 and 10" 3 L o . 



most of the ambiguous sources in GCMM probably lie at 
the near heliocentric distance. It is to notice that the last 
points toward lower luminosities could be considered to 
be the only few sources at the far distance, what would 
bring them within the range of luminosities assumed in 
the model (10~ 8 < L < 10- 3 L Q ). 

It is in principle possible to eliminate the near-far 
distance problem by selecting appropriate subsamples of 
sources from GCMM. We have then two possibilities: ei- 
ther we select the sources in the outer Galaxy (90° < I < 
270°) or we select all sources at the tangent point of the 
molecular ring (±50° < I < ±20°). The luminosity dis- 
tributions resulting from these selections were not better 
defined than the ones shown in Fig. [7] This is mainly due 
to the fact that the number of sources we are left with 
after the selection (m 90) is severely reduced as compared 
to the total. Also, the selection of sources at the tangent 
points does not completely eliminate the near-far ambi- 
guity, because of the difficulty of unambiguously defining 
the tangent point. 

It is important to notice that our estimates of the index 
of the luminosity function of methanol masers (—1.5 > 
a > —2) is in slight disagreement with what found in 
previous studies (~ —2). The difference in index could 
probably be reduced assuming a more complex luminosity 
function, as e.g. a broken power-law. We still consider the 
luminosity function presented here to be the best estimate 
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with the least number of assumptions and free parameters, 
and, most importantly, obtained without any selection of 
sources. The definitive shape of the luminosity function 
of methanol masers in the Galaxy will be found once the 
MMB Survey is completed. 

4. Conclusions 

Using the General Catalogue of 6.7 GHz Methanol Masers 
in the Milky Way, we have modelled their spatial distri- 
bution in the Galaxy and estimated their luminosity func- 
tion. We conclude the following: 

— methanol masers are distributed in a ring of some 5 kpc 
in radius around the Galactic centre; 

— most of the sources in GCMM are probably confined 
in a ring of heliocentric radius equal to i?o or slightly 
larger. This also means that the sources showing he- 
liocentric distance ambiguity lie probably at the near 
distance; 

— the luminosity function of methanol masers in our 
Galaxy is modelled as a power-law between sharp cut- 
offs having an index a between —1.5 and —2; 

— given a uniform survey which is able to define a 
methanol maser luminosity function, the approach pre- 
sented here will allow us to estimate the total number 
of sources in the Galaxy. As an example, a large scale 
survey with a sensitivity limit of 0.5 Jy will be able to 
detect w 50% of the total number of sources or more, if 
the power of the luminosity function is —1.5 or lower. 
This is what is expected from the MMB Survey. 
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